Two-dimensional periodic frustrated Ising models in a transverse field 
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We investigate the interplay of classical degeneracy and quantum dynamics in a range of periodic 
frustrated transverse field Ising systems at zero temperature. We find that such dynamics can 
lead to unusual ordered phases and phase transitions, or to a quantum spin liquid (cooperative 
paramagnetic) phase as in the triangular and kagome lattice antiferromagnets, respectively. For the 
latter, we further predict passage to a bond-ordered phase followed by a critical phase as the field is 
tilted. These systems also provide exact realizations of quantum dimer models introduced in studies 
of high temperature superconductivity. 
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The introduction of quantum fluctuations to a frus- 
trated system with a large ground-state degeneracy can 
lead to a multitude of new quantum phases and transi- 
tions. This occurs because the set of degenerate states 
provides no energy scale and hence all perturbations are 
singular. Any linear combination of the classically degen- 
erate states is a candidate quantum ground state. The 
most celebrated instance of such behaviour is the quan- 
tum Hall problem, where the degeneracy induced by a 
frustrated kinetic energy leads to a plethora of strongly 
correlated states when lifted by interactions and disorder. 

In this Letter we report results from the first system- 
atic study of a series of simple spin models which have 
the potential for such behaviour. These are geometrically 
frustrated transverse field Ising models. The frustrated 
nature of the exchanges leads to an extensive zero-point 
entropy, <S, in the classical (zero-field) case (for a review, 
see Ref. [[!)), and they are not ordered even at T = 0. 
The quantum dynamics is provided by the application of 
a magnetic field transverse to the couplings. 

These models are interesting as they are the simplest 
possible frustrated quantum magnets, containing the 
minimal ingredients required for a non-trivial quantum 
dynamics; this simplicity allows considerable progress 
towards their solution. Further, in addition to the en- 
couraging likelihood of direct experimental realization in 
highly anisotropic magnetic systems, we expect that they 
will serve as effective theories in other frustrated systems 
where a local Ising degree of freedom can be identified. 
Indeed, there is intense current interest in these models 
in the context of Heisenberg systems: the square lattice 
fully-frustrated Ising model arises in a group of theories 
of the cuprate superconductors g] . 

Specifically, we study Hamiltonians of the form H : 
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where the Jy are the nearest-neighbor exchange cou- 
plings With \Jij\ = J and riplaquette( J u7 J ) = "1. T is 

the strength of the transverse field, the S a are the Pauli 



spin operators and h is the strength of a (classical) lon- 
gitudinal, Ising symmetry breaking field. 

Counterintuitively, but commonly, there is an Ising 
version of "order by disorder" a la Villain B, in which 
configurations with long-range order allow the softest 
fluctuations and hence minimize the entropy/zero point 
energy; in the Ising case it is fluctuations into the ground 
state manifold that matter. In the complementary "dis- 
order by disorder" scenario, suggested by Fazekas and 
Anderson [pi, a disordered classical manifold continues 
into a disordered quantum phase. To our knowledge, no 
such disorder-free Ising spin liquid was known to exist 
previous to the work described here. 

To make progress in this strongly interacting problem, 
we have developed or adapted a number of approaches 
- variational, weak- and strong-coupling ones as well as 
mappings to other models. From our systematic study, 
we present three models realising different connections 
between classical and quantum ordering and, among oth- 
ers, unconventional critical and spin liquid phases. A 
more complete account of this study, including a wider 
range of lattices and finite-temperature properties, will 
follow this Letter M . 

Among our results, we derive a sufficient condition for 
quantum ordering, which is fulfilled by a class of mod- 
els, critical at T = 0, including the triangular lattice 
Ising (IAFM) antiferromagnet and the square fully frus- 
trated Ising magnet (FFIM). The kagome IAFM and the 
hexagonal fully frustrated Ising magnets, both exhibit- 
ing exponential correlations in the ground state manifold 
[Gj, do not fulfill this criterion and yet lead to differ- 
ent quantum states: the latter orders thus ruling out 
a strict connection between classical and quantum dis- 
order. By contrast, the former is a rare example of a 
quantum spin liquid in d = 2 + 1, as it remains disor- 
dered at all couplings. In the kagome problem there is 
a particularly interesting twist in which a purely longi- 
tudinal field, < \h\ < 4J, can be used to produce a 
critical classical finite-entropy phase which then orders 
in a further transverse field. Interesting byproducts of 





FIG. 1. Left: Maximally fiippable spin configuration of 
the triangular IAFM. Dimers mark the frustrated bonds, cov- 
ering the dual lattice. Flipping the spin in the centre, which 
experiences zero net exchange field, yields the indicated dimer 
plaquette move. Right: The FFIM on the hexagonal lattice, 
which can be thought of as a rectangular lattice (with vectors 
indicated by arrows) with a basis of four (encircled) sites. 
Broken lines represent antiferromagnetic bonds, solid lines 
ferromagnetic ones. The dimer configuration for the maxi- 
mally fiippable state with all spins aligned is shown. The 
elementary dimer move corresponds to flipping the two spins 
indicated by crosses. 



the analyses include unusual phase transitions, such as 
an 0(4) transition for the hexagonal FFIM. 
General considerations: Before discussing specific lat- 
tices, we sketch three general approaches. At small T, 
there is an Ising analog of the saddle-point method to 
identify instances of order by disorder. In this limit we 
need to diagonalize the transverse field term within the 
subspace of the classical ground states. Note that this 
will lead to a discontinuous change in the ground-state 
properties, even though the energy levels are set by T, 
and hence evolve continuously. As the quantum ground 
state is not perturbatively constructible, it is useful to 
think variationally about different configurations. Evi- 
dently, spins able to orient in the x-direction can gain 
energy from the transverse field. Such spins are "fiip- 
pable" in the z— representation, being connected to their 
neighbors by equal numbers of satisfied and unsatisfied 
bonds (see Fig. [y) . An intuitive variational candidate for 
the ground state is thus constructed from the classical 
ground states which maximize the number of fiippable 
spins; the privileging of these configurations, which are 
typically more regular, is the quantum mechanism of or- 
der by disorder. While this kind of semiclassical argu- 
ment is suggestive, it leaves open the possibility that the 
true wavefunction has the bulk of its support elsewhere 
in configuration space. Note that in cither case there are 
fluctuations in the ground state. 

A very useful feature of the transverse field problem is 
that at large T, there is a unique paramagnetic ground 
state (polarized along x) separated by a gap A = T from 
the lowest excited states, independent of wavevector. 
This degeneracy is lifted in an expansion in J/T: to first 
order, there appears a dispersion e(k) ~ A + ^\ e ik '-> Joj, 
which softens with reduced coupling. Even without high- 
order computations, one can identify soft wavevectors. 



Following Ref. 0, these can be used together with lat- 
tice symmetry considerations to guess at the Landau- 
Ginzburg- Wilson (LGW) action governing the transition 
and the symmetry-breaking pattern of the putative weak- 
coupling ordered phase. This possibility does not always 
exist; e.g. an XY magnet is a non-trivial problem on its 
own, which frustrates attempts to perturb in the Ising 
exchange. 

Finally, we derive an ordering criterion based on the 
mapping of the <i-dimcnsional transverse field Ising model 
onto a ferromagnetically stacked model in d + 1. This 
follows from the Euclidean path integral representation, 
Z = Tre~ H , of the problem governed by the 3D clas- 
sical Hamiltonian, H = J2{ij), n K tj^i( na r)S Z j{na T ) + 
^2 i n K T S*{na T )Sf((n+\)a T ), where a T is the imaginary 
time step introduced to obtain a discrete representation, 
Kf, oc Jij and the K T > are ferromagnetic. Time 
continuum quantum evolution corresponds to the double 
limit, a T — > 0, K T — ► oo with a T e KT /2 = T held fixed 
For small K T , the specific free energy, F, of this model 
is given by -0F = S + K*/2 E^(0)S,(0)) 2 + 0(K T 4 ), 
with \K S \ = oo enforcing the ground-state constraint. 

If the integral I = J d 2 r{SoS r ) 2 diverges, the free en- 
ergy above is non-analytic as K T — > 0, implying that the 
quantum (K T — ► oo) and classical points (K T — 0) are 
in different phases. For classically disordered antiferro- 
magnets with exponentially decaying spin correlations at 
T = even though no single term in the series expansion 
becomes unbounded, the whole series may diverge when 
K T is beyond the radius of convergence, and hence both 
ordered and disordered quantum states are allowed. In 
the following we demonstrate that all three possiblities 
are realised in practice by presenting an example for each 
one. 

Triangular IAFM: Order from criticality. Because 
this lattice is bond-sharing, the ground states of the tri- 
angular IAFM can be mapped onto hard-core dimer cov- 
erings of its dual hexagonal lattice (Fig. m; furthermore, 
since the latter is bipartite, the ground-state manifold of 
this spin model admits a height representation, which, 
under reasonable assumptions, implies that the classi- 
cal correlations are critical ||. We have a diverging 
Itri = J ol 2 r/r [|lj. Our above argument thence implies a 
transition separating the classical and the quantum mod- 
els. 

More directly, we study the effect of small T (h = 0) 
for this problem by noting that, in dimer language, 
states connected by the action of the resulting Hamil- 
tonian (Eq. m differ by a 60° rotation of a triplet of 
dimers (Fig. Ty). The degenerate problem is exactly the 
hexagonal-lattice version of the Rokshar-Kivelson quan- 
tum dimer model |l 

H QDM - -t (|V)(A| + h.c) + v (|V)(V| + |S)(A|) 

at the point t — T and v — pQ |. The ground-state 
of Hqdm at v = t is known to be an equal-amplitude 



superposition of all dimer coverings O]. As the di- 
agonal equal-time correlations at this point are exactly 
those of the classical problem, the question of fluctuation- 
selection reduces to determining if the point v = lies in 
the same phase. From an exact diagonalization study of 
clusters of this extended problem with up to 18 spins ||, 
we find that the the point v = lies in a hexagonal dimer 
phase which terminates at v = t. This is also the phase 
identified by our heuristic semiclassical analysis (Fig. nf). 

Such hexagonal correlations were previously found by 
an LGW analysis of ferromagnetically-stacked triangu- 
lar lattices J7|, which in the proper limit can be carried 
over to our model fl2|l . This analysis suggests that the 
paramagnetic phase orders via a 3D XY transition. 

We note that the presence of an XY-transition should 
manifest itself in a finite-temperature critical phase 
where T acts as a locking field to suppress defect pro- 
duction ||. Finally, we remark that these results apply, 
mutatis mutandis, to the square lattice FFIM. 
Kagome IAFM: Disorder from disorder. The 
Kagome IAFM has short-ranged correlations in the clas- 
sical limit, with a ground state entropy per site of 
0.502fcs, in excess of the triangular lattice value of 
0.323fcs IJJ. Its ground states map onto neither a height 
nor a dimer model. 

Interestingly, a longitudinal field \h\ < A J does not 
fully eliminate this extensive entropy but reduces it to 
precisely the triangular lattice result even as it creates a 
net ferromagnetic moment. Most remarkably, it simulta- 
neously induces critical (connected) spin correlators. To 
see this, note that maximizing the polarization requires 
that every elementary triangle of the lattice contain two 
up and one down spin. By placing dimers on the dual 
lattice defined by the triangles such that one dimer is 
centered on each down spin (Fig. |2j), one obtains dimer 
coverings of a hexagonal lattice precisely as in the trian- 
gular lattice problem in zero field, from which the above 
results follow. 

Next, consider the effect of tilting the field somewhat, 
i.e. of adding a small transverse component r -C h. At 
the first five orders in T/h there is no mixing between 
the polarized states and no lifting of their degeneracy. At 
sixth order they mix, and the degenerate manifold Hamil- 
tonian is precisely given by Hqdm with t oc T(T/h) 5 and 
v = 0. Hence we find that long-range bond order of 
the hexagonal kind considered previously sets in imme- 
diately. 

Consider now the problem with the transverse field 
alone. A low-order LGW analysis as described above fails 
to predict an ordering pattern: an entire dispersionless 
branch of excitations [pi goes soft, rather than only exci- 
tations at a small number of wavevectors. Following our 
previous strategy we can alternatively attempt to identify 
candidate orderings by finding the maximally flippable 
states. These states are precisely the 'dimer states' se- 
lected by a longitudinal field (and their Ising reversed 





disordered 



FIG. 2. Phase diagram for the kagome Ising antiferro- 
magnet in a field. Inset: Mapping of the kagome IAFM in 
a longitudinal field onto the hexagonal lattice dimer model. 
The down spins are marked by dimers. The non-trivial move 
to sixth order in perturbation theory corresponds to flipping 
all the spins marked by crosses. 



counterparts) since only they saturate the upper bound 
of two flippable spins in each triangle. A significant re- 
striction to this dimer manifold should manifest itself in 
correlations reminiscent of the bond-ordered phase. 

To test this possibility we have carried out Monte Carlo 
(MC) simulations on the Euclidean 3D path integral rep- 
resentation. As we expect any order to be strongest at 
small r, we have primarily simulated the problem with 
the spatial couplings made very large, in effect requir- 
ing that all time slices remain in the ground state man- 
ifold. In order to explore large K T without freezing, we 
use a cluster algorithm jbl in the time direction so that 
it is the exponential growth of the correlation length in 
the time direction that limits our simulations in practice. 
The largest system we simulated is a stack of height 512 
and the maximal coupling was chosen such that a typical 
cluster size is around 50. Each slice contains 432 spins. 

Fig. [3| shows the spin-spin correlation function of this 
system, which is evidently tiny beyond the first few neigh- 
bors, in stark contrast to the models known to order; e.g. 
in the triangular case, the saturated correlation function 
reaches well above 0.5 at the largest distances. 

Note that the shape of the quantum correlation func- 
tion mirrors that of the classical one: their signs are the 
same at all points in the region plotted. Since the clas- 
sical correlations are known to decay exponentially, and 
since the ratio of the two also appears to do so but with 
a smaller slope (inset, Fig [|), we conclude that the quan- 
tum correlations, although enhanced over the classical 
ones, still decay exponentially to zero with distance. 

In addition, we have checked that relaxing the spatial 
couplings leads to a further weakening of the correla- 
tions of the stacked model. We have also simulated the 
effect of a weak longitudinal field added to the transverse 
field and find that, while it induces a magnetization per 
site maft, the connected correlations appear to remain 




FIG. 3. Angularly averaged correlation functions (clas- 
sical and quantum) for the kagome lattice from MC simula- 
tions. Inset: The logarithm of the absolute value of the same 
functions. In addition, the points denote the logarithm of the 
ratio of the classical and the quantum correlation functions. 
The straight line is a guide to the eye. 



ing the presence of an ordering transition as T is lowered 
p5| . The connection of this ordering pattern with the 
variational analysis is rather involved ||; here we note 
that the fairly non-trivial ordering pattern is remarkable 
both as an instance of order by disorder and of ordering in 
a quantum dimer model where the classical correlations 
are short-ranged. 

In summary, we have found a rich variety of behavior 
in frustrated transverse field Ising models. The exciting 
possibility here is that it should be possible to approx- 
imate the triangular and kagome IAFMs fairly well not 
only in real, highly anisotropic magnets but also in more 
general quantum frustrated systems. 

We are grateful to J. Chalker for suggesting work on 
the kagome magnet. We thank D. Husc for several stim- 
ulating discussions. We would also like to acknowledge 
conversations with C.L. Henley and S. Sachdev. RM and 
SLS were supported in part by grants from the Deutsche 
Forschungsgemeinschaft, NSF grant No. DMR-9978074, 
the A. P. Sloan Foundation and the David and Lucille 
Packard Foundation. 



short-ranged. At short distances, we observe incipient 
correlations characteristic of the hexagonal dimer state 
described above for L <C h. This leads us to conjec- 
ture that, as the field is tilted further, there is likely a 
continuous quantum phase transition to the hexagonal 
phase. For larger h, diverging equilibration times pre- 
vent us from exploring this portion of the phase diagram 
numerically The simplest phase diagram incorporating 
all these facts within the framework discussed here is de- 
picted in Fig. 0; evidently, this is a fruitful topic for future 
work. 

Hexagonal lattice FFIM: Order from disorder. As 
an example of another classically disordered model, we 
consider the FFIM on the hexagonal lattice. This model, 
which maps onto a dimer but not a height model, is ob- 
tained from a ferromagnet by changing the sign of the 
exchange interaction of one bond in each hexagon, as 
displayed in Fig. Q. As the coordination of the lattice is 
odd, the transverse field does not lift the degeneracy un- 
til second order in L; this corresponds to flipping not an 
individual spin but a pair of neighboring spins (Fig. Py). 
Whereas the entropy of the ground states is extensive, at 
0.214fcs per spin H, the number of maximally flippable 
states (for an example, see Fig. m), is exponential only in 
the linear size of the system. 

We have carried out a LGW analysis as described 
above, and found that there are four soft modes, at 
wavevectors ±(71-/6, 7r/2), ±(5tt/6, 7r/2), with a surpris- 
ingly large unit cell of 48 spins. The resulting LGW- 
Hamiltonian |q] is 0(4)-symmetric up to sixth order, 
where a (in d = 2 + 1, dangerously irrelevant) 48- fold 
symmetry breaking term appears. 

Our MC simulations of the stacked model clearly dis- 
play the pattern predicted by the LGW analysis, indicat- 



R. Liebmann, Statistical Mechanics of Periodic Frus- 
trated Ising Systems (Springer, Berlin, 1986) and refer- 
ences therein. 

S. Sachdev and M. V ojta, hond-mat/991023^ , T. Senthil 
and M. P. A. Fisher, \:ond-mat/991022\ 
For work on one-dimensional models, see B.K. Chakra- 
bati et al., Quantum Ising Phases and Transitions in 



[2] 
[3] 

[4] 

[5] 
[6] 
[7] 

[8] 
[9] 

[10] 

[11] 
[12] 



[13] 

[14] 
[15] 



Tranverse Ising Models, (Springer- Verlag, Berlin, 1996). 
J. Villain et al, J. Phys. - Paris 41, 1263 (1980); E. F. 
Shender, Sov. Phys. JETP 56, 178 (1982). 
P. Fazekas and P. W. Anderson, Phil. Mag. 30, 23 (1974). 
R. Moessner, S. Sondhi and P. Chandra, in preparation. 
D. Blankschtein et al, Phys. Rev. B 29, 5250 (1984). 
B. Nienhuis et al, J. Phys. A 17, 3559 (1984). 
D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 
2376 (1998). 

C.L. Henley in A. Gervois, M. Gingold and D. Iagolnitzer, 
eds. Book of Abstracts of STATPHYS 20 (IUPAP Com- 
mission on Statistical Physics, Paris, 1998). 
M.E. Fisher and J. Stephenson, Phys. Rev. 132, 1411 
(1963). 

Alternatively, the quantum problem can be mapped onto 
a 3D height model; there, only smooth surfaces are ex- 
pected, which in turn suggests the observed ordering. We 
are grateful to S. Sachdev for explaining this to us. 
J. N. Reimers et al, Phys. Rev. B 43, 865 (1991) 
R. H. Swendsen et al, Top. Appl. Phys. 71, 75 (1992) 
This ordering pattern is already present for only 30 stacks 
of 1152 spins each. While our numerics are highly sug- 
gestive, a definitive finite-size scaling study is out of the 
reach of our current workstation due to the large size of 
the unit cell. 



